Density matrix renormalization group approach for many-body open quantum systems 
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The density matrix renormalization group (DMRG) approach is extended to complex-symmetric 
density matrices characteristic of many-body open quantum systems. Within the continuum shell 
model, we investigate the interplay between many-body configuration interaction and coupling to 
open channels. It is shown that the DMRG technique applied to broad resonances in the unbound 
neutron-rich nucleus 7 He provides a highly accurate treatment of the coupling to the non-resonant 
scattering continuum. 
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The DMRG method was introduced to overcome the 
limitations of Wilson-type renormalization group to de- 
scribe strongly correlated ID lattice systems with short 
range interactions (see recent reviews Exten- 
sions of DMRG to 2D lattice models are in general less 
accurate and require a reduction to a ID problem by 
selecting a path in the 2D plane. More recently, by re- 
formulating the DMRG in a single-particle (s.p.) basis, 
several applications to finite Fermi systems like molecules 
0, superconducting grains quantum dots Q, and 
atomic nuclei Q have been reported. The main advan- 
tage of DMRG is that it can accurately treat configu- 
ration spaces of dimensions well beyond the limits of 
a large-scale diagonalization. The early application of 
DMRG in the context of the nuclear shell model (SM) in 
the M-scheme exhibits some convergence problems that 
are not yet understood Q. The first DMRG+SM appli- 
cations in the angular-momentum-conserving J-scheme 
were successfully carried out in Ref. || and also re- 
cently in Ref. [lfj. A similar formalism, known as non- 
Abelian DMRG, has been previously applied to quantum 
lattice models 11.]. While most of the DMRG studies 
were focused on equilibrium properties in strongly cor- 
related closed quantum systems with hermitian density 
matrix, non-equilibrium systems involving non-hermitian 
and non-symmetric density matrices can also be treated 
with a DMRG procedure I n this work, we shall 

present the first detailed test of the DMRG approach to 
many-body open quantum systems (OQS) described by 
the Gamow Shell Model (GSM) [IlBj, i.e., the con- 
tinuum SM in the complex fc-plane. The DMRG ap- 
proach within the GSM is characterized by complex- 
symmetric reduced density matrices. We are particularly 
interested in developing an efficient strategy for selecting 
most important configurations involving s.p. resonances 
and non-resonant scattering states that describe many- 
body unbound states. Our approach was briefly sketched 



in Ref. |9|. 

GSM is the multi-configurational SM with a s.p. ba- 
sis given by the Berggren ensemble |l5j consisting of 
Gamow (resonant or Siegert) states and the non-resonant 
continuum of scattering states. The resonant states 
are the generalized eigenstates of the time-independent 
Schrodinger equation, which are regular at the origin 
and satisfy purely outgoing boundary conditions. GSM 
provides a quasi-stationary formalism for describing the 
time-dependent processes such as the mult i- particle de- 
cays. The s.p. Berggren basis is generated by a finite- 
depth potential, and the many-body states can be ex- 
panded in Slater determinants spanned by resonant and 
non-resonant s.p. basis states 0,0]- The GSM, which 
has been up to now applied to weakly bound/unbound 
atomic nuclei, can be equally applied for the description 
of other OQS, e.g., open microwave resonators and self- 
bound atomic systems, such as 3 He,/v neutral droplets at 
the limits of their stability 0|. Another possible appli- 
cation of the GSM are quantum dots, where the interplay 
between electron-electron correlations and the continuum 
coupling yields the transition in the conductance proper- 
ties. 

At present, the principal limitation of GSM applica- 
tions is the explosive growth in the number of config- 
urations (i.e., dimension of the many-body Fock space) 
with both the number of active particles and the size of 
the s.p. space. To ensure completeness of the Berggren 
basis, for each resonant s.p. state of quantum num- 
bers lj, with I being the orbital angular momentum 
and j the total angular momentum, one should include 
a large set of discrete non-resonant continuum states 
{(Zj-)W; £ = 1, . . . , M} c l ying on a continuous contour L l ^_ 
in the complex fc-plane [13l Il4j . These continuum states 
become new active shells in the many-body framework of 
GSM and, because of their presence, the dimension of the 
GSM Hamiltonian matrix grows extremely fast. This ex- 
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plosion of the Hilbert space is much more severe than in 
the standard SM, which deals solely with a limited num- 
ber of discrete s.p. states. The use of the Berggren en- 
semble implies a complex- symmetric Hamiltonian matrix 
for the hermitian Hamilton operator. Due to the large 
number of non-resonant continuum shells whose proper- 
ties vary smoothly along the contour, the GSM Hamil- 
tonian matrix is significantly denser than that of a con- 
ventional SM. As a result, matrix operations are more 
time-consuming in GSM than in SM, and the Lanczos 
procedure for GSM eigenvectors 0] becomes very slow. 

Most of the configurations involving many particles in 
the non-resonant continuum contribute very little to low- 
energy GSM eigenfunctions calling for a smart selection 
of the most important configurations. For that purpose, 
we propose a new algorithm based on the DMRG proce- 
dure for finding the GSM eigenvalues. The main idea is to 
gradually consider different s.p. shells of the discretized 
non-resonant continuum in the configuration space and 
retain only N opt optimal states dictated by the eigen- 
values of the density matrix with largest modulus. The 
procedure will be illustrated using the example of many- 
body resonances in the neutron-unbound nucleus 7 Hc de- 
scribed in terms of N v =3 active (valence) neutrons out- 
side the closed core of 4 He. 

The s.p. basis is generated by a Woods-Saxon (WS) 
potential with the radius Ro=2 fm, the depth of the cen- 
tral potential Vo=47 MeV, the diffuseness <i=0.65 fm, 
and the spin-orbit strength V so =7.5 MeV 0. This 
s.p. potential reproduces experimental energies and 
widths of the s.p. resonances 3/2]" (g.s.) and 1/27/ 
(first excited state) in 5 He at £7=0.745-i0.32MeV and 
£7=2.13 — i2.94MeV, respectively. The neutron valence 
space consists of the Opi/2 and 0p 3 / 2 resonant states and 
the corresponding non-resonant states {pi/2}c, {P3/2}c- 
The L P + /2 contour in the complex-fc plane (in the fol- 
lowing, k is expressed in units of fm -1 ) is defined by 
a triangle with vertices at (TZe(k),lm(k))= (0,0), (0.33,- 
0.33), (0.5,0.0), and a segment along the 7£e(fc)-axis from 
(0.5,0) to (1.0,0.0). Similarly, L+ 3/2 is given by a triangle: 
(0,0), (0.17,-0.17), (0.5,0), and a straight segment from 
(0.5,0) to (1.0,0). Each segment of these contours is dis- 
cretized with the same number of points corresponding to 
the abscissas for Gauss-Legendre quadrature. The GSM 
Hamiltonian is a sum of the WS potential, representing 
the effect of an inert 4 He core, and the two-body interac- 
tion among valence neutrons. The latter is approximated 
by a finite-range surface Gaussian interaction [l7^ with 
the range (J.— 1 fm and the coupling constants depending 
on the total angular momentum J of the neutron pair: 
V (J=0) =-542 MeV fm 3 , V (J=2) = -479 MeV fm 3 . These 
constants are fitted to reproduce the binding energies of 
6 He and 7 He with respect to the core. 

For the purpose of DMRG, the configuration space 
is divided into two Fock subspaces: A (built from the 



s.p. resonant shells 0pi/2, 0p 3 / 2 ), and B (built from 
s.p. non-resonant shells {p\/2~}c and {j>3/2}c)- In the 
warm-up phase, one calculates and stores all possible ma- 
trix elements of suboperators of the Hamiltonian in A: 
at, (at a) K , (a^a^) K , ({a^a^) K a) L ), {(a^ a^) K (aa) K )), and 
constructs all the states \k)A with n = 0, 1, • • • ,N V par- 
ticles coupled to all possible j-values. Then, from the 
sets: {pi/2}c, {pi/2}c one picks up the first pair of s.p. 
states, calculates the full set of matrix elements of subop- 
erators for this added pair, and constructs all the states 
\i) with n particles (n = 0, 1, ••• ,N V ). The subspace 
B is gradually enlarged by successively adding pairs of 
non-resonant shells (pi/2 , P3/2) until the number of states 
\i)b exceeds iV pt- Then the Hamiltonian is diagonalized 
in the space {\k) A \i) B } J made of vectors in A and B. 
Obviously, the total number of particles in such states 
must be kept equal to the total number of valence par- 
ticles, and the angular momentum J is conserved. To 
find the many-body resonance from all the eigenstates, 
\ty J ) = ^2cki{\k) A \i) B } J , one selects the eigenstate hav- 
ing the largest overlap with the GSM eigenstate in a pole 
approximation (i.e., diagonalized in the subspace A only). 
From this eigenstate, one calculates the reduced density 
matrix p u } 3 = J^k c fci c fci' with a fixed value of j in states 
\i) b, 10 b - By construction, the density matrix p B is 
then block-diagonal in j. Within the metric defining the 
Berggren ensemble, the density matrix is complex sym- 
metric; hence, its eigenvalues are complex. The reduced 
density matrix is diagonalized, and N op t eigenstates \u v ) 
having the eigenvalues with the largest absolute values 
are retained. Those eigenvalues correspond to the most 
important states in the enlarged set. All the matrix ele- 
ments of suboperators for the optimized states are recal- 
culated. Then, the next pair of non-resonant continuum 
shells is added and, again, only N ovt eigenvalues of the 
density matrix are kept. This procedure is repeated until 
the last shell in B is reached, providing a first guess for 
the wave function of the system in the whole ensemble of 
shells forming the basis. From now on, all the s.p. states 
are considered. 

At this point, the warm-up phase ends and the so- 
called sweeping phase begins. Here, one constructs states 
with 0, 1, 2, • • • , N v particles and then the process contin- 
ues in the reverse direction until the number of vectors 
becomes larger than N opt . If the m th pair of p-shells in 
B is reached, the Hamiltonian is diagonalized in the set 
of vectors: {\k,i prev )\i)} J , where i pre v is a previously op- 
timized state (constructed from first m — 1 pairs of shells 
in £?), and i is a new state. The density matrix is then 
diagonalized and the N opt i-states are kept. The pro- 
cedure continues by adding the (m — l) st pair of shells, 
etc., until the first pair of p-shells in B is reached. Then 
the procedure is reversed again: the first pair of shells is 
added, then the second, the third, etc. 

In the examples presented in this work, the number of 
s.p. shells included in A and B are N T =2 and N c = 60 
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FIG. 1: GSM+DMRG results for the real (top) and imaginary 
(bottom) parts of the g.s. (J n = 3/2~) energy of 7 He plotted 
as a function of the number of DMRG steps. The number 
of states kept during the sweeping phase is A r opt =20 (dotted 
line), 22 (dashed line), and 24 (solid line). 



(30 shells in each continuum: {pi/2j-c, {P3/2jc), respec- 
tively. In the warm-up phase, 8 vectors are kept. Fig- 
ure n shows the convergence properties of the ground 
state (g.s.) J 71 ' = 3/2~ resonance energy in 7 He plot- 
ted as a function of the number of DMRG steps A^tep- 
(As discussed earlier, one step corresponds to adding a 
pair of {pi/2} c , fe3/2jc shells.) Converged results, man- 
ifested by a local plateau at the complex energy E=Eq, 
are found after N s t cp < 30 steps. Then, with increasing 
Agtcp, one enters the unstable region (40 < N < 70). 
A second plateau, corresponding to practically the same 
energy Eq as the first plateau (and fairly independent of 
Aopt), appears at A/step ~ 90. The "plateau- unstable" 
sequence repeats periodically. It is worth noting that the 
deviation from TZe(Eo) can be both positive and negative. 
This is because the Hamiltonian matrix is not hermitian 
but complex-symmetric; hence, the usual reasoning based 
on the Ritz variational principle does not apply. 

The behavior of the DMRG procedure, as can be seen 
in Fig.^ an d Tabled shows a gradual reduction of the en- 
ergy variation around the plateau with increasing number 
of sweeps N sw . Most importantly, it decreases quickly 
with the number of states -/V opt kept in the sweeping 
phase. The relative precision of the DMRG procedure 
is extremely high: ~ 10~ 6 for A r opt =24 and A r sw =4. 

Figurc[2]illustrates the convergence of the DMRG pro- 
cedure for the first excited state in 7 He, which is calcu- 
lated to be a broad resonance with T=3.37 MeV. Here, 7 
vectors are kept in the warm-up phase. The full conver- 
gence is attained in about 2 sweeps, i.e., somewhat slower 
than for the narrow g.s. resonance (r=0.14 MeV). Still, 
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TABLE I: Relative precision of the real part of the g.s. energy 
in 7 He as a function of the number of states N op t kept in the 
sweeping phase at the second (A r sw =2) and fourth (A r sw =4) 
sweep of GSM+DMRG. 
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FIG. 2: Similar to Fig. Q except for the first excited state 
(J n = 1/2") in 7 He. 



the resulting plateau is excellent. Again, as in Fig. ^ 
i?o is practically independent of -Aopt. These two exam- 
ples demonstrate that the proposed DMRG algorithm is 
extremely efficient in selecting the most important non- 
resonant continuum structures. 

A major problem in studies of OQS is to maintain the 
completeness of the many-body basis containing contri- 
butions from discretized non-resonant s.p. continua. In 
this respect, the Berggren ensemble is not different: it 
has to be discretized for any practical applications of 
the GSM. The accuracy of discretized many-body cal- 
culations relics upon the independence of calculated ob- 
servables such as energies, transition probabilities, and 
cross sections from the number of s.p. shells in the com- 
plex non-resonant continuum. To illustrate this point, we 
show in Fig. |3the g.s. energy of 7 He as a function of the 
total number of s.p. shells N s h- One can see that the de- 
pendence oilZe(E) and Xm(E) on A/ S h is not monotonic, 
showing oscillations which disappear only if a sufficiently 
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FIG. 3: Convergence of the real (top) and imaginary (bottom) 
parts of the g.s. energy of 7 He as a function of the the total 
number of shells N sh (= N c + N T = N c + 2) included in the 
s.p. basis. As in Fig. Q 8 vectors are kept in the warm-up 
phase. In the sweeping phase, 7V pt=22 states are kept. The 
inset shows the imaginary part of the g.s. energy of 6 He as a 
functiuon of iV s h. Here, N opt —4 (10) vectors were kept in the 
warm-up (sweeping) phase. 



dense discretization is applied. A similar oscillatory pat- 
tern is seen in the imaginary part of the g.s energy of 6 He 
(shown in the inset). Since 6 He is bound, the deviation of 
Tm(E) from zero is due to too coarse discretization; the 
completeness is practically achieved for iV s h=70. While 
the error due to an insufficient number of non-resonant 
shells is of the order of 1CU 2 keV, the accuracy of DMRG 
is in the range of 10~ 6 keV. 

The results shown in Figs. Ill.il demonstrate that 
the fully converged GSM results with respect to both 
the number of sweeps and the number of shells in 
the discretized continua can be obtained using the 
GSM+DMRG algorithm. The rank of the biggest ma- 
trix to be diagonalized in GSM+DMRG grows extremely 
slowly with N^. In the example shown in Fig. [3J it 
changes from d=941 to <i=1001, whereas the dimension 
of the GSM matrix in the J-coupling scheme varies as 
D oc AT s 3 h , from L>=6149 for iV sh =26, to D=332,171 for 
Af sh =98. Hence, the gain factor D(N c )/d(N c ) ~ iV s 3 h 
quickly grows with increasing N c (or the size of block 
B). 

Our results indicate that (i) the complex energy Eq 
at the DMRG plateau very weakly depends on AT opt , 
and (ii) the precision of calculated complex energies 
weakly depends on N s h, provided that N s h is sufficiently 
large. This means that the convergence features of the 
GSM+DMRG procedure can be tested by varying iV s h 



and A^opt. Once those parameters are optimized in 'small- 
scale' GSM+DMRG calculations, the final calculations 
can be performed in the large model space to obtain fully 
converged results. 

In summary, this study describes the first application 
of the DMRG method to unbound many-fcrmion systems 
described by a complex-symmetric eigenvalue problem. 
The encouraging features of the proposed GSM+DMRG 
approach open the scope to systematic and high-precision 
studies of the continuum coupling effects in multiparticle 
OQS such as open microwave resonators, multi-electron 
open quantum dots, atomic nuclei, and atomic clusters 
close to the particle drip line. 
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